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ABSTRACT 

We present unpublished Spitzer IRAC observations of the HH 1/2 young stehar outflow processed 
with a high angular resolution deconvolution algorithm , that produces sub-arcsecond (~ 0.6 - 0.8'^) 
images. In the resulting mid-infrared images, the optically invisible counterjet is detected for the flrst 
time. The counterjet is approximately half as bright as the jet at 4.5/im (the IRAC band that best 
traces young stellar outflows) and has a length of ~ 10'^ The NW optical jet itself can be followed back 
in the mid-IR to the position of the exciting VLA 1 source. An analysis of the IRAC colors indicates 
that the jet/counterjet emission is dominated by collisionally excited H2 pure rotational lines arising 
from a medium with a neutral Hydrogen gas density of ~ 1000 — 2000 cm~^ and a temperature ~ 1500 
K. The observed jet/counterjet brightness asymmetry is consistent with an intrinsically symmetric 
outflow with extinction from a dense, circumstellar structure of ~ &^ size (along the outflow axis), 
and with a mean visual extinction, Ay~ 11 mag. 

Subject headings: circumstellar matter — stars: formation — ISM: jets and outflows — infrared: ISM 
— Herbig-Haro objects — ISM: individual objects (HHl/2) 



1. INTRODUCTION 

HH 1 and 2 were the flrst detected HH objects (Herbig 
1951; Haro 1952). Since then, they have played an im- 
portant role in the study of outflows from young stars, 
particularly because most of the general characteristics 
of HH outflows were flrst seen in HH 1 and 2 (see the 
review of Raga et al. 2011). 

The flrst measurement of proper motions in HH objects 
(Herbig & Jones 1981) showed that HH 1 and 2 formed 
part of a bipolar outflow. The source of this bipolar out- 
flow (centered between HH 1 and 2) was discovered in the 
radio continuum by Pravdo et al. (1985). This "VLA 1" 
source was later shown to have a jet-like structure (of 
~ V^ X 0'^2, see Rodriguez et al. 2000), aligned with the 
HH 1/2 axis. 

A jet-like structure directed towards HHl (the "HHl- 
jet") is observed at optical wavelengths (Strom et al. 
1985 point out this feature, which was also visible in older 
images of the region). The base of this jet-like structure 
approaches the position of the outflow source (VLA 1) 
in images taken at progressively longer IR wavelengths 
(Roth et al. 1989). 

Reipurth et al. (2000) presented optical and IR images 
of the HHl-jet region obtained with the HST. Their NIC- 
MOS H2 2.12/im and [Fe II] 1.64/im images show that 
the emission of the HHl-jet extends to within ~ 2^^ from 
the VLA 1 outflow source. It is well known that all these 
multiwavelength jets are a manifestation of the same as- 
trophysical phenomena (see e.g. the review by Raga et 
al. 2010). The fact that the observed HHl-jet emission 
does not extend to the position of the VLA 1 source and 
that a counterjet (directed towards HH 2) is not detected 
appears to be due to the presence of a dense molecular 
structure approximately centered on VLA 1 (see Tor- 
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relies et al. 1994; Choi & Lee 1998 and Cernicharo et al. 
2000). 

In this paper we present new Spitzer observations in 
which the counterjet (directed from the VLA 1 source 
towards HH 2) is detected for the flrst time, and ending 
nearly thirty years of speculation about the nature of its 
absence. 

2. OBSERVATIONS 

The observations of HH 1/2 system come from the 
GTO Spitzer Space Telescope program by Giovani Fazio 
(PID 43) on the 'Orion Streamers' obtained with the 
infrared camera IRAC (Fazio et al. 2004) in February 
2004. The data, consisting of the basic calibrated frames 
or BCDs, have been recovered from the Spitzer Legacy 
Archive, version S8.18. The original surveyed area is 
large and covers approximately 0.75°x3.3°in the four 
bands (1, 2, 3, 4) = (3.6, 4.5, 5.8, 8.0) /im. The data were 
collected using the High-Dynamic-Range (HDR) mode 
with a 12sec integration time for the 'long' frames and 
0.6sec for the 'short' ones. For this work we have used 
the 'long' frames only, with a total integration time of 
48sec per pixel. 

The BCDs were then reprocessed with the HiREs de- 
convolution software AWAIC (A WISE Astronomical Im- 
age Co-Adder) developed by the Wide Field Infrared 
Survey Explorer (WISE) for the creation of their Atlas 
images (see e.g. Masci & Fowler 2009) P] The AWAIC 
software optimizes the coaddition of individual frames 
by making use of the Point Response Function (PRF) 
as an interpolation kernel, to avoid flux losses in under- 
sampled arrays like those of IRAC, and also allows a 
resolution enhancement (HiRes) of the flnal image, by re- 
moving its effect from the data in the deconvolution pro- 
cess. A similar method has been applied on the Spitzer 
data of young stellar outflows like HH 46/47 (Noriega- 
Crespo et al. 2004a; Velusamy et al. 2007) and Cep E 
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(Moro-Martm et al. 2001; Noriega- Crespo et al. 2004b; 
Velusamy et al. 2011), quite successfully. On IRAC im- 
ages, the HiRes processing enhances the angular resolu- 
tion from the standard - ^2!' to - 0.6'' - 0.8'' (Velusamy 
et al. 2007). 

In Figures 1 and 2 we show the images of the HH 1/2 
system obtained in the four IRAC bands (3.6, 4.5, 5.8 and 
8.0/im) before (i.e. standard coaddition of the BCDs) 
and after the HiRes reprocessing (iteration 25), respec- 
tively. At the iteration 25 the HiRes AWAIC process- 
ing reaches an optimal angular resolution, simultaneously 
preserving most of the structure of the surrounding dif- 
fuse emission, for instance that of the arc near the South 
of the images. Some small artifacts, however, are present 
in all bands, indicating that not further iterations are 
needed. Figure 3 shows a closer look of the 4.5/im im- 
age, where the newly discovered counterjet is the bright- 
est, marking as well some of the well known optical knots 
of the HH 1 and 2 objects, and the position of the VLA 1 
source. The bright infrared source along the outflow sym- 
metry axis, the so called Cohen- Schwartz stars (C-S), 
that once was thought to be the driving source (Cohen 
& Schwartz 1979) is also indicated, plus the VLA 2 source 
that drives the HH 144 outflow (Reipurth et al. 1993) 

One can appreciate the performance of the HiRes pro- 
cessing and the ability of the mid-IR observations to dis- 
cover new features by comparing the processed IRAC 
4.5/im with archival data at optical and near-IR obser- 
vations of the HH 1/2 system obtained by the Hubble 
Space Telescope^ with ~ 5 times better angular resolu- 
tion (Fig. 4). The optical image was taken by the Wide 
Field Planetary Camera 2 (WFPC2) in 2007 (Hartigan 
et al. 2011) at [S II] 6717/31 A, while the near-IR was 
taken with the NICM0S3 camera in 1998 at v=l-0 H2 
2.12/im emission line (Reipurth et al. 2000). It is quite 
remakable how similar are the vibrational (at 2.12/im) 
and rotational H2 emission that arises from the S(9), 
S(10), S(ll), S(12) transitions at 4.952, 4.408, 4.180 and 
3.996/im, covered by IRAC 4.5/im band. About halfway 
between the VLA 1 source and HH 2, there are two or 
three condensations detectable at 2.12 and 4.5/im, and 
given their positions could belong to counter flow as well. 

3. THE HH 1 JET AND COUNTERJET 

A closer look of the [SII] 6717/31 A, 2.12, 4.5 and 
8.0/im emission around the jet/counter region is shown 
in Fig. 5. A detailed analysis of the optical and near- 
IR HH 1-jet properties has been carried by Reipurth et 
al. (2000), where they showed how the H2 2.12/im and 
[Fe II] 1.64/im emission arise closer to the VLA 1 source, 
only 2.5" NW of it, than that of [S II]. Their analysis on 
the extinction, using the fact that the [Fe II] 1.54/im and 
[S II] 0.67/imhave similar excitation energies and ioniza- 
tion potential, and that their ratio is nearly constant for 
weak shocks, allowed them to determine a 4 mag increase 
in the ^ 5" that the nir-IR emission gets closer to the 
VLA 1 source. 

In the mid-IR at 4.5/im both the jet and counterjet 
are clearly detected, and both can be traced back to the 
VLA 1 source. The jet /counterjet path, unfortunately, 
is very close to two bright infrared sources; o n the NE, 
the C-S source and at ^ 10" SE of VLA 1 by a fainter 
one. And so it is possible that very close to the circular 




Figure 1. IRAC maps of the HH 1/2 system, from top left clock- 
wise, at 3.6, 4.5, 5.8 and 8.0/im; using a inverse grayscale where 
dark regions represent high intensity. The field-of-view is ~ 5^ 
North is up and East is left. 
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Figure 2. IRAC observations of the HH 1/2 system as in Fig 1, 
after 25 iterations using the HiRes AWAIC algorithm. 



edges were the AWAIC HiRes algorithm suppresses the 
sources, that couple of knots are affected by this artifact. 
Conservatively, the jet has a ~ 12.9^' length in the NW 
direction arising from VLA 1, while the counterjet has a 
~ 10^' length in the opposite direction. Both the optical 
and the mid-IR jet extend out to the same NW position, 
the A^- knot (Eisloffel, Mundt & Bohm 1994). Anf fi- 
nally, at both 4.5 and 8.0/im couple of arcseconds NW 
of VLA 1, the emission broadens following the structure 
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Figure 3. A closer view of the H 1/2 system at 4.5/im, after the 
HiRes AWAIC processing, where the jet and counterjet are best 
detected. The names of some of the weh known optical knots of 
HH 2 (SE) and HH 1 (NW) are included (see e.g. Raga, Barnes &; 
Mateo 1990), as well as those of HH 144 flow (Reipurth et al. 2000) 
its driving source (VLA 2), and the bright IR Cohen-Schwartz (C- 
S) source (Cohen Sz Schwartz 1979). 

observed at the base of the 2.12/im jet, the so called 'X' 
nebula (Reipurth et al. 2000). This nebula seems to be 
connected with the knots East of the jet; they do have a 
similar structure at 2.12 and 4.5/im. 

3.1. The Medium Surrounding the Counterjet 

In Figure 6, we present 4.5 and 5.8/im intensity trac- 
ings along the jet /counterjet system. In order to obtain 
these tracings, we have defined an axis parallel to the di- 
rection of the outfiow (the x-axis of Figure 6, with x = 
corresponding to the position of the VLA 1 source, and 
positive X towards the NW), and averaged the intensity 
in a direction perpendicular to the outfiow, in a box ex- 
tending ±3'^ to each side of the outfiow axis. A back- 
ground was computed from contiguous boxes on each 
side of the outfiow (of l'^5 widths), and subtracted at 
all positions x from the axial box in order to obtain the 
jet emission. 

The intensity tracings show that the emission from 
the NW jet is generally stronger than the counterjet 
emission. The bottom frame of Figure 6 shows the 
jet /counterjet intensity ratios as a function of distance 
from the VLA 1 source at 4.5 (left; broken line) and 5.8 
(right; solid line). It is clear that the jet/counterjet in- 
tensity ratio is > 1 at 5.8/im(at all positions). At 4.5/im, 
the jet/counterjet intensity ratio is > 1 out to ~ 8^' from 
the source, and < 1 at ~ 10^' from VLA 1. This re- 
gion (with jet/counterjet intensity ratio < 1) is associ- 
ated with the last knot along the counterjet, seen close to 
the stellar source ~ 12^' to the SE of VLA 1 (see Figures 
3 and 5), and might therefore be contaminated by the de- 
convolved point spread function of the star. Therefore, 
we conclude that the observations are consistent with 
jet/counterjet intensity ratios > 1 at all positions both 
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Figure 4. A comparison of the HH 1/2 system, from left to right 
at optical (WFPC2 [SH] 6717+31 A), near-IR (NICMOS 2.12/im) 
and mid-IR wavelengths (IRAC 4.5/im HiRes [25 iterations] pro- 
cessing). The position of VLA 1 source is marked with a ?>" in 
diameter circle. 




Figure 5. A closer view of the jet/counterjet region, using once 
again the Hubble Space Telescope observations in [SH] and 2.12/im, 
and compared with those at IRAC 4.5 HiRes (25 iterations) pro- 
cessing. The position of VLA 1 source is marked with a 1" in 
diameter circle. 



at 4.5 and 5.8/im (at least out to ~ 8^' from the VLA 1 
source, see Figure 6). If one assumes that the jet and 
counterjet emission is intrinsically symmetric, one would 
attribute the position-dependence of the jet/counterjet 
intensity ratio (/j//cj, see the bottom frame of Figure 
6) to the extinction produced by a high density clump 
surrounding VLA 1. The fact that Ij/Icj reaches a peak 
(with a value of ^ 4.5 at 4.5/im, and of ?^ 3 at 5.8/im, see 
Figure 6) at x ~ 3'^ would then indicate that the high 
density clump has a projected diameter of ~ 6'^ and 
that for larger values of x the counterjet emerges from 
behind the clump (lowering the observed Ij/Icj ratio). 
This ~ &' size is consistent with the extension along the 
outflow axis of the flattened H^^CO+ clump observed by 
Choi & Lee (1998), surrounding VLA 1. 

The formation of stellar outflows arising from young 
stellar objects (YSOs) are expected to be symmetric, 
since there is a preferential rotational axis and to the best 
of our understanding of the formation of proto-stellar jets 
requires a magnetohydrodynamical coupling of the ac- 
creted gas with the spherical symmetric protostar, likely 
to have also a symmetric bipolar magnetic structure (see 
e.g. Pudritz et al. 2007, and references therein). The ac- 
creted gas is provided by a disk-like structure created as 
a result of the original cloud spinning during the proto- 



stellar collapse (see e.g. Klein et al. 2007, and references 
therein). One could imagine, however, that the accretion 
process does not need to be symmetric, and could feed 
the protostar in a time dependent alternative way (e.g. 
one star pole at the time). If this was the case, then 
the symmetry of the stellar jet could be broken, at least 
over a certain period of time, with no material ejected in 
one direction or another. And if this was the case then 
one could trace back the steps to understand where (e.g. 
given the morphology of the protostar's magnetic field, or 
its coupling with the accretion disk or the transfer of gas 
and angular momentum, etc) this symmetry is broken. 
A true asymmetric young stellar jet could have profound 
consequences in our understanding of the low mass star 
formation process. For nearly thirty years we have won- 
dered if the HH 1/2 system had a counterjet, although 
as mentioned in the introduction evidence for very dense 
gas (n(H2 = 10^ cm~^) structure around VLA 1, using 
Ammonia as a tracer (Torrelles et al. 1994), suggested 
that extinction was playing a major in hidding that com- 
ponent of the flow. 

The NH3 observations have an angular resolution of 
~ 4" (~ 2 X 10^ AU at the distance of Orion), sampling 
very well a 2' x3' region covering the HH 1/2 outflow, 
and enough to distinguish a pancake-like structure per- 
pendicular to the outflow axis and a East-West tempera- 
ture gradient in it that indicates a further asymmetry. It 
is no a simple disk-like structure, but certainly a dense 
structure consistent with the idea of "a collapsing in- 
terstellar (~ 0.4pc) toroid around VLA 1" (Torrelles et 
al. 1994). NH3 has been used as tracer of high density 
and low temperature molecular gas around young stellar 
outflows for more than two decades. The NHS (1,1) tran- 
sition requires H2 densities higher than 5xl0^cm~^ to 
be detected reliably (Torrelles et al. 1983). Typical val- 
ues around YSOs range from n(H2) = 5xl0^cm~^ to 
~ lO^cm"^ and Temperatures of 15 - 30 K (Torrelles et 
al. 1983). The structure around the VLA 1 in HH 1/2 
has n(H2) ~ 10^ cm~^and 12 K. Similar conclusion ws 
reached by Choi & Zhou (1997) using other high density 
molecular tracer (HCO+). Slightly better angular reso- 
lution observations (4.3^^ x 2. 8'^ beam, Choi & Lee 1998) 
have traced even closer to the VLA 1 source the toroidal 
structure. Very close to the VLA 1 source, within a 
diameter of 4 AU, Cernicharo et al. 2000 estimated 
a very high visual extinction (80-100 magnitudes), but 
they were let's concerned with the "toroidal structure". 
In Figure 7 we show a comparison of the NH3(1,1) obser- 
vations by Torrelles et al. (1994) with the mid-IR emis- 
sion at 4.5/im. Two things are worth noticing. First, that 
high density molecular gas does indeed prevent us from 
detecting the counterjet at shorter wavelengths, and sec- 
ond, that the morphology of this dense structure is likely 
to prevent us to detect the HH 144 counterflow. 

Independently of detailed structure of the "toroid", it 
is clear that a large column density of material is a long 
the line-of-sight of the counterjet, and this is higher than 
in line-of-sight of the visible jet. An estimate of the ex- 
tinction around the counterjet is possible, although re- 
quires some necessary assumptions. For instance, based 
on symmetry one could assume that the emission prop- 
erties at 3.6 and 4.5/mum of jet and counterjet are very 
similar over the same area, and that any difference be- 
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Figure 6. Top frame: surface brightness (averaged across the 
width of the jet, in MJy/sr) at 4.5 (dashed hne) and 5.8/im (sohd 
hne) along the NW jet (positive x) and SE counterjet (negative x) 
centered on the VLA 1 source, showing the asymmetry of the two 
outflow lobes. Bottom frame: the jet/counterjet intensity ratio as 
a function of distance from the VLA 1 source at 4.5 (left panel; 
broken line) and 5.8 /im (right panel; solid line). 

tween the jet (North) and the counterjet (South) are due 
only to extinction. A comparison of the ratio of the sur- 
face brightness of the 3.6 to 4.5/mum images, indeed 
shows smaller values on the South with respect to the 
North. The average ratio of the North and South sec- 
tions of the jet, using a 2.4''x39.0'' box, give values of 
0.40±0.27 and 0.20±0.16, respectively, i.e. the counterjet 
is twice as faint when compared to the jet. Two recent 
studies have analyzed the properties of the mid-IR ex- 
tinction (Indebetouw et al. 2005, Cambresy et al. 2011), 
in regions located in the Galactic plane, i.e. RCW49 and 
the Trifid, with very similar results, so one can assume 
that such extinction approximately holds in the environ- 
ment around the HH 1/2 system. Taking the ratios of 
extinction at Visual, K5 (2.2/im), 3.5 and 4.5/im to be 
Ax/Ay = 0.112, As.e/Ax = 0.611 and A4.5/AK = 0.500 
(Cambresy et al. 2011), and neglecting the extinction at 
4.5/mum (since both jet/counterjet are detected) then at 
3.6/im at the counterjet one obtains, A3.6 = 0.75, that 
corresponds to A^ = 1.2 and Ay = 11.0 magnitudes, 
respectively. 

4. THE HH 1/2 MID-IR KNOTS 

Since some of the first IRAC images of stellar outflows 
from Spitzer (e.g. Noriega-Crespo et al. 2004a, b, Morris 
et al. 2004), it has been clear that the bands at 4.5 and 
5.8/im are particularly suitable to study them. There 
are seven pure rotational H2 emission lines, from S(12) 
at 3.996/im to S(6) 6.107/imthat fall within their band- 
passes. The 8/im band includes some bright v=0-0 H2 
lines (e.g. S(5) 6.907/im), but the emission in this band 
tends to be dominated by broad band features from very 
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Figure 8. IRAC color-color diagram (after Ybarra &; Lada 2009) 
of compact knots in the HH 1/2 outflow with optical counterparts, 
including HH IF (IF), HH 2 (A, E, F, G, H, I, K, and Z) and 
regions around the jet (Jet), counterjet (CtrJet) and the VLA 1 
source (see Table 1). 

of bowshock-like condensations. For example, knots like 
HH IF, HH 2A, HH 2G and HH 2H are high excitation, 
while HH 21, HH 2L are low excitation (Raga, Bohm, & 
Canto 1996). For the time variable jet, where the relative 
velocity between knots determines essentially the shock 
velocity, the overall excitation is known to be low and 
the shocks to be weak (Reipurth et al. 2000). Stronger 
shocks will produce higher post-shock gas densities and 
temperatures, but could also dissociate H2; stronger H2 
emission indeed is expected to arise from low velocity 
shocks (either J-type or C-type (see e.g. Draine & Mc- 
Kee 1993, for a review). 

So how the properties of the mid-IR knots, dominated 
by emission of H2 rotational lines at 4.5 and 5.8/im, com- 
pare with those of the optical knots? One of the highest 
optical excitation knots, HH IF, is slightly above the 
4000 K model, the region of the diagram where indeed 
H2 dissociates. The jet and counterjet, low excitation 
regions, do fall within the low temperature ~ 1500K and 
low density (~ 10^ cm~^) models. Notice also that their 
IRAC colors are very close to each other, providing sup- 
port to the assumption that their emission properties 
are nearly identical. The low excitation HH 2F and the 
region around VLA 1, are closer to the low 1500 K tem- 
perature model, but at higher gas density ~ 10^ and 
10^ cm~^, respectively. The rest of the mid-IR knots 
are found between these two extremes. We notice that 
HH 2H and G, that are along the highest density model 
(10^cm~^) are also colder than the low excitation knots 
HH 21 and L, and this suggests that faster cooling at 
higher densities could play a role in reducing the excita- 
tion in the mid-IR. 

5. SUMMARY AND CONCLUSIONS 

We present previously unpublished Spitzer observa- 
tions of the HH 1/2 outflow. The outflow is detected in 
the four IRAC channels (3.6, 4.5, 5.8 and 8.0/im). In the 
4.5/im channel, a well collimated jet /counterjet system 
is seen emanating from the VLA 1 radio source, which 
represents the discovery of the SE counterjet (previous 
optical and IR observations only detected the NW jet 

^ http://irsa.ipac. caltech.edu/data/SPITZER/docs/irac/iracinstrumei^j^^^^o|<Q^Q^2:ds HHl. 



Figure 7. A comparison of the dense structure detected in 
NH3(1,1) at 10.5 km s~i (dark contours; from Torrelles et al. 1994) 
with the emission at 4.5/im around the HH 1/2 jet /counterjet re- 
gion. 



small UV stochastically heated dust particles or Poly- 
cyclic Aromatic Hydrocarbons (PAHs) (see e.g. Tie- 
lens 2008). Recently, Ybarra & Lada (2009) have used 
a combination of IRAC colors to uniquely identify the 
thermal emission arising from collisionally excited H2 as 
a function of the HI gas density and temperature, and 
thus provide a reliable way to study the thermal struc- 
ture of the stellar outflows by using the ([3.6] — [4.5]) and 
([4.5] -[5.8]) colors. 

Some of the knots of the HH 1/2 system are clearly 
compact and at a first approximation one can treat them 
as point sources and transform their photometric fluxes 
into magnitudes (see e.g. "IRAC Instrument Handbook" 
F] , §4.11.1, 'Best Practices for Extended Sources'). One 
needs to be aware that interpreting surface brightness 
measurements into colors can be off by 5% - 10% ("IRAC 
Instrument Handbook", §4.11.1). We adopt a 10% un- 
certainty across bands in our measurements, and use the 
standard IRAC zero points ("IRAC Instrument Hand- 
book", Table 4.1) of F^o = 280.9, 179.7, 115.0 and 64.9 
Jy for the 3.6. 4.5, 5.8 and 8.0/im bands, respectively. 
The results are presented in Table 1 and display in a 
color-color diagram in Figure 8, following Ybarra & Lada 
(2009). The diagram shows the properties of shocked ex- 
cited H2 as a function of constant HI density (dotted line) 
and temperature (solid line). At a temperature higher 
than 4000 K, the H2 is likely to be dissociated (Ybarra 
& Lada 2009). 

The excitation properties of the HH 1/2 knots at op- 
tical and UV wavelengths are very well known (see e.g. 
Solf, Bohm & Raga 1998; Solf et al. 1991; Bohm & Solf 
1992; Schwartz et al. 1993; Bohm, Noriega-Crespo & 
Solf 1993; Moro-Martin et al. 1996; Molinari & Noriega- 
Crespo 2002) and are the result of their shock velocity 
and relative geometry. In general, high excitation knots 
are compact and interact through strong shocks (60 - 100 
km s~^) with the surrounding medium, while low excita- 
tion emission arises from weaker shocks or the "wings" 



We find that the ratio between the jet and counter- 
jet emission is strongly dependent on distance from the 
VLA 1 source. The jet/counterjet ratio is ~ 2 times 
larger for the 4.5/im than for the 5.8/im emission. This 
result is consistent an extinction effect (which would 
be stronger at shorter wavelengths). The fact that the 
jet/counterjet ratio peaks at a distance of ~ 3^' from the 
VLA 1 source indicates that the extinction is probably 
produced by a dense structure having an extent of ~ &^ 
along the outflow axis. Such a dense molecular struc- 
ture around VLA 1 has been observed by Choi & Lee 
(1998) and Cernicharo et al. (2000). We also show that 
the large structure of the dense gas traced by NH3. plays 
role in hidding both the HH 1/2 and HH 144 counterjets. 
Under the assumption that the ratio of 3.6 to 4.5/im is 
constant for the jet/counterjet, we estimate a Ay = 11 
mag. 

Finally, by using the IRAC colors ([3.6] - [4.5]) and 
([4.5] - [5.8]) and assuming that outflow emission is dom- 
inated by shocked excited H2 (Ybarra & lada 2009), we 
show that some of the compact mid-IR knots share sim- 
ilar excitation properties as those of determined for the 
optical knots (Raga, Bohm & Canto 1996). 
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Table 1 

HH 1/2 System IRAC Colors 



Knot 


[3.6] - [4.5] 


[4.5] - [5.8] 


RA 


Dec 


AperRad(^0 


Z^ 


1.29±0.13 


1.51±0.16 


5h36m27.40s 


-6d47m22.2s 


3.0 


F 


1.73± 0.18 


0.58± 0.06 


5h36m26.84s 


-6d47ml5.0s 


1.8 


L 


0.70± 0.07 


1.72± 0.17 


5h36m26.65s 


-6d47m28.3s 


3.0 


E 


1.48± 0.15 


1.15± 0.12 


5h36m26.58s 


-6d47ml3.2s 


1.8 


K 


1.36± 0.14 


1.11± 0.11 


5h36m26.24s 


-6d47m06.1s 


1.8 


G 


1.25± 0.13 


0.93± 0.10 


5h36m26.00s 


-6d47ml0.1s 


1.8 


A 


1.04± 0.10 


1.53± 0.15 


5h36m25.34s 


-6d47mll.9s 


2.4 


H 


0.93± 0.10 


0.92± 0.10 


5h36m25.65s 


-6d47ml7.2s 


2.4 


I 


0.79± 0.08 


1.26± 0.13 


5h36m25.53s 


-6d47m22.5s 


1.8 


Jet 


2.46± 0.25 


1.00± 0.11 


5h36m22.44s 


-6d45m58.1s 


1.8 


CtrJet 


2.40± 0.24 


1.10± 0.10 


5h36m23.00s 


-6d46m08.0s 


1.8 


VLA 1 


2.04± 0.20 


0.84± 0.09 


5h36m22.83s 


-6d46m05.0s 


1.8 


IF 


0.71± 0.07 


2.03± 0.20 


5h36m20.23s 


-6d45m05.6s 


2.4 



Knot 8 at 2.12/im (Davis, Eisloffel &; Ray 1994) 



